Structures and Dielectric Properties of Amorphous High-K Oxides: Hf02, Zr02 and their alloys 
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High-K metal oxides are a class of materials playing an increasingly important role in modern device physics 
and technology. Here we report theoretical investigations of the properties of structural and lattice dielectric 
constants of bulk amorphous metal oxides by a combined approach of classical molecular dynamics (MD) - 
for structure evolution, and quantum mechanical first principles density function theoiy (DFT) - for electronic 
structure analysis. Using classical MD based on the Bom-Mayer-Buckingham potential function within a melt 
and quench scheme, amorphous structures of high- k metal oxides Hfi_a;Zra;02 with different values of the con- 
centration X, are generated. The coordination numbers and the radial distribution functions of the structures are 
in good agreement with the corresponding experimental data. We then calculate the lattice dielectric constants 
of the materials from quantum mechanical first principles, and the values averaged over an ensemble of samples 
agree well with the available experimental data, and are very close to the dielectric constants of their cubic form. 

PACS numbers: 61.43.Er, 77.55.D-, 77.84.-s 
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I. INTRODUCTION 

In the complementary metal oxide semiconductor (CMOS) 
devices, Si02 has been the first and only choice of gate ox- 
ide material for a long time. However, with the continued 
decrease in the feature size of the CMOS devices, Si02 is 
no longer reliable as a gate oxide due to the high tunnel- 
ing leakage current through it at small thickness. Major re- 
search efforts are continuously being carried out in search for 
a suitable replacement of the Si02 as a gate materialji Sev- 
eral promising candidates are metal oxides such as Hf02f^ 
Zr02!^ — and Al203j^ all of which have high value of dielectric 
constant k. High-K metal oxides in their amorphous (a-) form 
are more preferable as a gate oxide over their crystalline form 
due to several important advantages they provide: (i) isotropic 
physical properties; (ii) no crystalline domain boundary which 
leads to less defects at the interface with the Si substrate; and 
(iii) good compatibility with the conventional CMOS fabri- 
cation process. Some alloy structures of these high-K metal 
oxides are also being studied extensivelji^^ and, in fact, a Hf 
based alloy material has already been in its third generation of 
production as a gate oxide in the semiconductor industry and 
further improvements of thermal stability, dielectric constant 
and material preparations are underway .^^iS 

In clear contrary to the abundance of experimental results 
in the literature, theoretical studies of the structure and di- 
electric properties of amorphous metal oxides and their al- 
loys are quite limited. One of the reasons is the difficulties 
in generating reasonable and reliable amorphous structures 
with the available theoretical methods. Experimentally, the 
high-K materials are deposited on silicon substrate by vapor 
deposition (as-deposited films) followed by annealing pro- 
cesses at around 1000 Kiiiii Simulating the deposition pro- 
cess and the resulting amorphous high-K material structures 
by molecular dynamics (MD) is extremely time consuming 
because the deposition rate must be controlled close to the ex- 
perimental value that is usually very low. The as-deposited 
a-Hf02 and a-A^Oa films have been generated by kinetic 
Monte Carlo (KMC) methods^iil with certain assumptions 



such as the atoms are kept fixed on pre-determined crystal 
sites and less important dynamics processes of atomic motion 
are omitted. Both MD and KMC are useful for the simulation 
of the initial atomic layer deposition processes of the amor- 
phous film growth. To generate bulk amorphous structures, 
ab inito MD simulations can be used in a melt and quench 
scheme, a-Hf02 r a-ZrO a"^'^'^^ and their silicates^ia have 
been successfully simulated this way. The activation-relation 
technique is another method to generate structures of contin- 
uous disordered systems^S and has been applied to generate 
a-Zr02 structuresJi In this method, one optimizes the struc- 
tures many times to determine a local energy minima (in con- 
figurational space) using a force field which can be calculated 
from quantum mechanical first principles or described by an 
empirical potential function. For the former the calculation 
cost is high while for the latter a reliable potential function is 
required. Recently, several extended-Tersoff potentials'^, were 
proposed for certain metal oxides and applied to generate their 
amorphous structures from classical MD 

Amorphous structures are long range disordered systems 
which typically requires much larger than ^ 100 atoms to 
simulate. Another issue, perhaps more difficult to deal with, is 
the simulation time scale that is needed for the atoms to evolve 
into the amorphous structure. While ab initio methods have 
simulated up to picoseconds time scale, one certainly wishes 
to investigate the structures and electronic properties of metal 
oxides at much larger time scales and much larger sizes in or- 
der to reveal the long time limit of the structural evolution and 
the physical properties of the resulting amorphous material. It 
is the purpose of this work to report such an investigation. 

In particular, we calculate the lattice dielectric constants of 
bulk amorphous metal oxides by a combined approach of clas- 
sical MD - for structure evolution, and quantum mechanical 
first principles DFT - for electronic structure analysis. We 
overcome the size and time scale difficulty by classical MD 
simulations in a melt and quench scheme and generate amor- 
phous structures of Hf02, Zr02 and their alloys in the form 
of Hfi_j;Zr2:02 with different values of the concentration x 
(0 < a; < 1). The classical MD is based on the interatomic 
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FIG. 1: (color online) Structures of a-Hf02 generated by classical 
MD starting from the cubic ciystal form of Hf02 with different lat- 
tice constants, corresponding to different densities. If the density 
is high, the system reciystallizes to its cubic form (a) or any other 
crystal form (b); if the density is low, the system tends to form disor- 
dered structure with big holes (d), whereas with moderate density of 
7.97 g/cm"^ which coixesponds to lattice constant of 5.6 A for crys- 
talline Hf02, a reasonable amorphous structure is generated (c). Blue 
spheres represent Hf atoms, red spheres are for O atoms. 



0.5 



O 0.4 






0.2 


'c 






0.15 






s 






0.1 






DF 


0.05 






0^ 



Coordination number 







Hf-Hf 










(b) 




HJ 















Distance (A) 



FIG. 2: (color online) (a) The coordination number and (b) the ra- 
dial distribution function (RDF) of the amorphous Hf02 structure 
obtained by the classical MD simulations. Green, red and blue lines 
in (b) represent Hf-O, O-O and Hf-Hf RDF, respectively. The results 
are obtained from structures having 768 atoms. 



Born-Mayer-Buckingham potential functioriiai This potential 
has been developed over many years and extensively applied 
to investigate the growth processes of Y2O3, Zr02 and Mg- 
Al-0 thin films.22i2S, Here we show that high-K bulk materi- 
als of amorphous nature can be generated quite efficiently us- 
ing this particular potential function. Detailed analysis of the 
amorphous structures are carried out including the coordina- 
tion number, the radial distribution function and the potential 
energy per atom of the structures. Using the atomic struc- 
tures generated this way, we carry out first principles density 
functional theory (DFT) calculations of the lattice dielectric 
constant and the results are averaged over ten independent 
samples for each material. The sample averaged dielectric 
constants show very good agreement with the corresponding 
experimental data. 

The combined approach of classical MD and DFT is an ex- 
cellent way to deal with the problem of large sizes and large 
time scales required by the materials physics of metal oxides. 
Several advantages can be summarized: (i) classical MD is in- 
expensive for simulating systems containing large number of 
atoms which is important for a proper characterization of the 
amorphous structures having long range disorders, here we 
have gone up to 768 atoms although several thousand atoms 
can be easily handled; (ii) in the melt and quench scheme a 
much longer relaxation time scale, here to several nano sec- 
onds, can be evolved after the material is heated above the 
melting point. Such a long time scale is important for eras- 
ing the memory of the initial structure thereby the final out- 
come becomes independent of the initial conditions. Here we 
show that the same kind of amorphous structure is obtained 
by evolving from different initial structures including from a 



completely random structure, (iii) Since the calculated physi- 
cal properties of amorphous material should be averaged over 
an ensemble of samples due to the randomness of the struc- 
ture, the combined classical-quantum approach allows such 
configurational average to be done. Here we show that the cal- 
culated dielectric constants of the metal oxides become spa- 
tially isotropic after configurational average, but it is typically 
anisotropic if only one sample is calculated. 

The paper is organized as follows. In Sec. II we present the 
classical MD method as employed in this work. In Sec. Ill we 
discuss the structural properties of our simulated amorphous 
structures, and also present our results for the lattice dielectric 
constants of these structures. Section IV concludes the paper 



II. METHODS 

In this section we discuss the potential function used in our 
classical MD simulations along with a short description of 
the melt and quench scheme. In our MD method the inter- 
actions between the atoms are described by the Born-Mayer- 
Buckingham potential function-' as shown below: 

V{r) = Aexp{-'^)~^ + -^. (1) 

p r°j iTTEonj 

The cation-cation interactions are assumed to be purely 
coulombic, the parameters for O^^ — O^^ interactions and 
cation-anion interactions are presented in Table. U All the pa- 
rameters are obtained from Ref. ll2l1l . We employed Wolf et. 
al. 's method^^'^** for the summation of the Coulomb interac- 
tion term. 
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At the initial stage of the simulation, all the atoms are as- 
signed velocities taken from a Maxwell's distribution at the 
room temperature T^. The system is then heated up to a tem- 
perature T,n = 4000 K, well above the experimental melt- 
ing point of around 3000 K, followed by an evolution of 5 
ns at the same temperature; and then quenched back to the 
room temperature Tq. During this melt and quench process, 
the motion of an atom is governed by the Langevin equation 
= F — av + Fr- Actually, the temperature of the sys- 
tem is adjusted by resetting the velocity of a randomly se- 
lected atom at the interval of Ai according to the following 
expression 
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where 6* is a parameter between and 1, v-^™ (^) is a veloc- 
ity selected randomly from the Maxwell's distribution at tem- 
perature Tm or Ta via a random number ^. Velocity resetting 
with Eq. Q corresponds to taking the damping coefficient 
a — m9/2lS.t and the random force Fr = mO^/'^v'^ / l\t 
in the Langevin equation, where m is the atomic mass.^^ Pre- 
vious calculations-"^showed that a value of 0. 1 for 9 is reason- 
able and appropriate. After the system is quenched to Ta, it 
is allowed to evolve for 10 ps to calculate the radial distribu- 
tion function (RDF). Finally, the system is cooled down to 
K by a damping method,'^ and thus concludes the procedure 
for generating one amorphous structure of a particular metal 
oxides. 

In the simulations, the initial structure is modeled with a 
2x2x2 super cell of Hfi_j;Zrj;02 containing 96 atoms; or 
4x4x4 containing 768 atoms. A periodic boundary condi- 
tion is imposed on the supercell. The initial atomic configu- 
ration is the cubic crystalline form or in the form of a random 
structure by randomly placing the constituents atoms (with the 
correct ratio) inside the simulation box. These drastically dif- 
ferent initial conditions verify that the MD simulation time 
scales given above are large enough since the same kind of 
final amorphous structures (similar dielectric constant values 
and potential energy per atom) are obtained. 

TABLE \: Parameters for Bom-Mayer-Buckingham potential func- 
tion as obtained from Ref. I2lll. 



elements charge A(eV) p(A) C 

Hf +4 1454.6 0.350 0~ 

Zr +4 1453.8 0.350 

O -2 22764.3 0.149 20.4 



III. RESULTS 

A. Structure properties of YiSi x'Ltx02 

To generate reliable amorphous structures it is important to 
find the optimal density to use for the initial structure. Hence, 
we first perform a series of independent simulations on the 
96-atom super cell, starting from initial structures of different 
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no. 3: (color online) Potential energies per atom for the amorphous 
structures of Hfi_iZra;02 averaged over 10 samples for each value 
of X (solid red line). The potential energies per atom for each indi- 
vidual sample are also shown in the figure (circles). The potential en- 
ergies of these samples reside in the range of —36.88 eV to —36.82 
eV, which are about 0.3 eV higher than that of the cubic Hf02. The 
very small differences in the potential energies among the samples 
with the same value of x confirms the reliability of the structures 
generated. 



lattice constants {i.e. different densities) of cubic Hf02 and 
Zr02 in the range from 4.5 A to 6.0 A with an interval of 0.1 
A. We found that structures of amorphous nature (see Fig.lc) 
are obtained when the lattice constant is set to 5.6 A which 
corresponds to the density of 7.97 g/cm'^ for Hf02 and 4.66 
g/cm-^ for Zr02. If we start with a lattice constant less than 
5.5 A, the systems recrystallize to a cubic form (Fig. la) or 
some other crystal form (Fig. lb); while with a lattice constant 
greater than 5.8 A, the systems tend to form disordered struc- 
tures with big holes inside (Fig. Id). Our optimized density 
values are close to those obtained previously from ab initio 
calculations which were 8.6 g/cm^ for Hf02^ and 4.86 g/cm^ 
for Zr02'^ Hence, the lattice constant of 5.6 A is used as 
initial condition (cubic form) in subsequent melt and quench 
evolution to the amorphous structures of Hfi_a;Zr2;02. 

To verify the generated amorphous structures we calculate 
the coordination numbers and the radial distribution functions. 
Results are presented in Fig.|2]for a-Hf02 obtained from the 
768-atom supercell. From Fig. we note that the three- 
coordinated O atoms and six-coordinated Hf atoms dominate 
the a-Hf02 structures which agree with the known ab initio 
results.ii^ From Fig. we find that the first peaks in the Hf- 
O, O-O and Hf-Hf RDFs are at 2.1 A, 2.8 A, and 3.6 A, re- 
spectively, again compare well with ab inito studiesi^aS These 
comparisons indicates that the classical MD can reliably gen- 
erate amorphous structures of the metal oxides. Importantly, 
the efficiency of classical MD allows us to generate ensem- 
bles of amorphous samples of Hfi_3;Zra;02 for each given 
concentration x, which is essential for statistical analysis of 
the dielectric properties (see below). 

For further analysis, we calculate the potential energy per 
atom of the a-Hfi_^Zrj;02 structures for different values of 
X ranging from to 1 with an interval of 0.125. For each 
value of X, we generated 10 different amorphous samples by 
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TABLE II: Calculated lattice dielectric constants for Hf02 and Zr02 
relative to the dielectric constant at vacuum. 
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FIG. 4: (color online) Lattice dielectric constants for the amor- 
phous structures of Hfi_a;Zr3;02. All the black symbols are for the 
values in each direction of X,Y and Z averaged over 10 samples, 
while the square symbols connected with a solid red line represent 
the values which are averaged over all the three directions for each 
value of a;. In the inset the standard deviations of the dielectric con- 
stant values over the 10 samples along the X direction have been 
presented. 



varying the initial conditions, and then calculated the potential 
energy per atom for each sample which are presented in Fig. |3] 
The results show that the potential energies of a-Hfi_2;Zr2;02 
structures as generated are located in the range from —36.88 
eV to —36.82 eV which is a reasonable and small spread - 
higher than that of their crystal forms, and have similar poten- 
tial energies among themselves. We also generate a structure 
of a-Hf02 starting from an initial system of 32 Hf atoms and 
64 O atoms placed randomly in a cubic simulation box. The 
calculated potential energy of this particular a-Hf02 struc- 
ture comes to —36.86 eV which is in the same energy range 
of the a-Hf02 structures that were generated from the cubic 
Hf02 initial conditions. It is very satisfactory that the gen- 
erated amorphous structures are not dependent on the initial 
condition, indicating that the simulation time scales presented 
above are adequate. 



B. Dielectric properties of Hf i _ a; Zr^ O2 

A most important property of the high-K metal oxides is 
their dielectric constant. In this section we present our calcu- 
lations of the lattice dielectric constant (static) of Hfi_a;Zrx02 
amorphous structures generated in the last section. We found 
that the electronic contribution to the dielectric constant is 
substantially smaller than the lattice contribution and does not 
change very much for different structures, hence it will not be 
discussed until when we calculate the total dielectric constant 
to compare with the measured data (see below). 

For the calculations of lattice dielectric constants we em- 
ploy density functional perturbation theory within a projector 
augmented wave method in the generalized gradient approxi- 
mation (GGA) using Perdew-Burke-Ernzerhof (PBE) param- 
eterization as implemented in the VASP software.-'^ The pseu- 
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dopotentials for Hf and Zr include only the outermost shells 
in the valence. An energy cutoff of 400 eV is used in the plane 
wave by sampling the Brillouin zone only at the T point. The 
MD-generated amorphous structures have not been relaxed 
with DFT in VASP. This is reasonable since we generated by 
MD a statistical ensemble of structures and the values of the 
dielectric constant have been averaged over the ensemble. For 
each amorphous structure the dielectric constant is calculated 
for an ensemble of 10 different samples (the 96-atom super 
cell). The sample averaged values are presented in Table 
for a-Hf02 and a-Zr02 along with the dielectric constants for 
all 3 of their crystal forms. We notice that the lattice dielec- 
tric constants of a-Hf02 and a-Zr02 are close to that of their 
cubic form, suggesting that some cubic-like short range order 
content may dominate the amorphous structures. Indeed, re- 
cently experimental measurement&ii on sputtering samples of 
a-Hf02 discovered a cubic-like short range coordination in the 
amorphous microstructure which led to high measured dielec- 
tric constant k ^ 30. Our calculated lattice k (see Table Hlb 
is 24.52, adding the electronic contribution which we found 
to be around 4.4, the total k w 28.9 which agrees quite well 
with the measured data. To confirm that the k of a-Hf02 being 
close to that of the cubic Hf02, is not due to the cubic initial 
condition of the structure simulation, we have evolved a single 
a-Hf02 structure from a completely random initial condition 
(see previous section): its lattice dielectric constant is found 
to be 21.785, 25.43 and 25.163 in the three spatial directions 
respectively; its electronic contribution is found to be 4.45. 
These values are totally consistent with that obtained from the 
cubic initial conditions. 

In comparison to previous theoretical first principles 
investigationS j V- our k values are higher This difference 
may be explained by noting that in Ref. IslfslflTll the amor- 
phous structures were obtained from a relaxed monoclinic 
structure with cubic lattice vectors which tends to have lower 
K values. We should mention that direct comparison with ex- 
perimental measurement is often difficult due to scarcity of 
data for dielectric constants of amorphous bulk structures . 

Note that the contribution of a given mode to the lattice 
dielectric constant scales as Z^'^/wf, where Z* is the mode 
effective charge and uj\ is the frequency of the Ath IR-active 
phonon mode, a large k will be obtained with the existence of 
modes with simultaneously high-Z* and low-o;.^'"* It is easy to 
understand that cubic Hf02 (Zr02) with modes of higher-Z* 
and lower-w has larger k value than monoclinic structure with 
modes of lower-Z* and higher-w. The amorphous structures 
may have the similar Z* and oj values with their cubic form. 
This consideration should provide a way to explore possible 
short range structures within the a-Hfi_a;Zra;02. 

The ensemble averaged lattice dielectric constants of all 
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the amorphous structures generated for different values of the 
concentration x are presented in Fig. |4] We observe a grad- 
ual increase (solid red line) in the values of dielectric constant 
from a-Hf02, through Hfi_j;Zrj;02, to a-Zr02. We also note 
that the spatial resolved dielectric constants in X, Y and Z di- 
rections (represented by the black symbols in Fig.|D are quite 
isotropic which is expected for amorphous structures. It is 
very important to note that an isotropic dielectric constant is 
obtained only when the values are averaged over the ensem- 
ble of samples of the same structure (same x). For a single 
sample, we always observe large spreads in the values of di- 
electric constant in different directions. This finding under- 
pins the major advantage of employing a rather inexpensive 
classical MD to efficiently generate amorphous structures for 
the statistical average of the physical properties. Note that all 
the values of dielectric constants presented in this study are 
relative to the dielectric constant at vacuum. 



IV. SUMMARY 

Using classical MD simulations within a melt and quench 
scheme, the amorphous structures of Hf02 and Zr02, along 
with their alloy structures of the form of Hfi_^Zr^02, have 



been generated. A suitable potential function, namely the 
Born-Mayer-Buckingham potential, has been adopted for 
these simulations. The calculated coordination numbers and 
radial distribution functions of the generated structures agree 
well with the previous studies. The values of the ensemble 
averaged lattice dielectric constants of these amorphous struc- 
tures are nearly isotropic, close to that of their cubic form, and 
compare well with available experimental data. On the other 
hand, a single sample always produces anisotropic k values. 
Our results indicate that the combined approach of classical 
MD - for structure evolution, and quantum mechanical first 
principles DFT - for electronic structure analysis, is a reliable 
technique for investigating high-K metal oxides. 
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